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Abstract 

The gravitational waves (GW) background generated in a double inflation- 
ary model, with two scalar fields mutually interacting through gravity only, 
and its relative contribution T / S to large-angle temperature fluctuations of the 
relic microwave background are investigated in detail. The relation between 
T/ S and the slope of the GW spectrum n-r is shown to be a discriminative test 
between a slow-roll inflation driven by one scalar field and more complicated 
models. It is found that the GW amplitude is not exactly zero in minima of 
spectral oscillations, this property being an observational, in principle, mani- 
festation of GW being in a squeezed vacuum state during inflation. 

PACS Numbers: 04.62. +v, 98.80. Cq 







Cosmological models in which the early Universe passed through two successive 
de Sitter (inflationary) phases of quasi-exponential expansion and other pa- 

pers) are usually called double inflationary ones. The initial spectrum of adiabatic 
perturbations in these models is not approximately fiat (i.e., n 1, or Harrison- 
Zeldovich like), it has a characteristic scale and more power on larger scales. This 
represents one of the possible ways to reconcile the cold dark matter model (CDM) 
with observations of large-scale spatial distributions of galaxies and their peculiar 
velocities ^ and large-angle fluctuations of the relic microwave background tempera- 
ture [0, because the observational data imply an excess of power in the perturbation 
spectrum on scales beginning from (25 — 50)h~^Mpc up to the present Hubble radius 
Rh = cHq^ = 3000/;.~^Mpc as compared to the CDM model with fiat initial spec- 
trum of adiabatic perturbations (the standard CDM model) normalized to produce 
a correct galaxy-galaxy correlation function and having the biasing parameter 6 ~ 2 
on a scale of 8/i~^Mpc. 

The most logically economic double inflationary models possess only 2 additional 
free dimensionless parameters as compared to the standard CDM model. There- 
fore, they still retain most of the predictive power of the inflationary scenario. De- 
tailed comparison of a number of double inflationary models with recent observational 
data |Tl|, |1^ shows that they are compatible with observations (though for a rather 
small window of parameters). In p|, we performed a detailed investigation of a sim- 
ple double inflationary model with two massive scalar fields and found initial spectra 
of adiabatic perturbations and gravitational waves generated from vacuum quantum 
fluctuations during a double inflationary stage analytically in the case when there is a 
pronounced intermediate stage between two stages of inflation. As explained in ||T2 



power spectra for which the relative height of the two "plateaus" in the spectrum 
of adiabatic perturbations is compatible with observations correspond to values of 



10 — 15, the latter number increases to 18 if the recently proposed fit to the 



COBE data corresponding to Qrms-ps = 20 for n = 1 ||T3| is used (this is only 
a little bit less than 15-25 that was proposed in p before the announcement of the 
COBE measurements). The intermediate stage between the two inflationary stages is 
not pronounced in that case. Therefore we have to resort to numerical simulations in 
order to get the perturbation spectra with good accuracy. A unique and remarkable 
property of some variants of the inflationary scenario (in particular the chaotic and 
power-law inflations) is that such a large primordial gravitational wave background 
is generated from vacuum quantum fluctuations that it is fairly reasonable that a 
modest, but still significant part of the large- angle fluctuations of the cosmic mi- 
crowave background (CMB) temperature ^'^^t' '^'^ observed by COBE |jlO[ is due to 



these gravitational waves. That is why it is very important to calculate the ratio of 
the gravitational waves amplitude to that of scalar perturbations and their relative 
contribution to the temperature anisotropy in the double inflationary scenario. 

In this letter, we extend this analytical and numerical study to the spectrum of 
primordial gravitational waves and their contribution to the large-angle temperature 
fluctuations of the CMB. Though the ratio of GW amplitude to that of adiabatic 
perturbations and the corresponding ratio of their contributions to the large-angle 



1 



AT/T fluctuations appear to be rather small in the case under consideration, this 
study gives a possibility to obtain important basic conclusions about a discrimina- 
tive observational test between slow-roll inflation driven by one scalar fleld and more 
complicated models, in particular double inflation, and about an observable (in prin- 
ciple) feature indicating that GW were in a pure squeezed state for some time during 
inflation. 

Let us give the essentials of our model of double inflation and generation of per- 
turbations in it. The phenomenological Lagrangian density describing matter and 
gravity during and near the inflationary stage, is given by 

L = + li<l^h,,K - mM) + liMf - m]<f,) (1) 

where = 0, .., 3, c = = 1, vrih ^ m; and the Landau-Lifshitz sign conventions are 
used. The space-time metric has the form 

ds^ = dt^ - a^{t)6ijdx'dx^, i,j = 1, 2, 3. (2) 

Spatial curvature may always be neglected because it becomes vanishingly small dur- 
ing the flrst period of inflation driven by the heavy scalar fleld. The background 
flelds are homogeneous and classical. The perturbed FRW background, in longitudi- 
nal gauge, has a metric given by 

ds^ = (1 + 2$)rft2 - a^{t){l - 2^)6ijdx'dx^ . (3) 

The same mechanism that produces quantum fluctuations of the scalar flelds will 
also produce tensor metric perturbations or gravitational waves. For each of the 
polarization states A, where A = 1,2, and the polarization tensor is normalized by 
the condition ej^e*-' = 2, the amplitude hx is given by 



hx = VlQnG(j)x (4) 
where (px is a massless scalar fleld satisfying 

□ = 0. (5) 
Let us consider the spectrum of tensor and scalar adiabatic perturbations 

h{k) = Bk-""^ , = j^Ak-^'^ ■ (6) 

The quantities A and B are those introduced in |15|, they are only weakly de- 
pending on k while h{k) and $(A;) are deflned as follows 

(Va, /^k'.A') = \h\k)5^^\\^ - k')5AA' , (7) 
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where the Fourier transform conventions are = — ^- f $fr)e '^^'^d^k. We assume 

also that a(t) oc t^^^ today. 

Let us first consider the ratio of the amphtudes of GW and adiabatic perturbations 
generated in double infiation. For wavelengths well outside the Hubble radius, | <^ if, 
each Fourier mode hk xij]) has a constant value. If they cross the Hubble radius during 
a quasi-de Sitter stage, h{k) is given by 

h{k) = Vl6^^^ (9) 



where tk is the Hubble radius crossing time defined as = H{tk). In (|^) both 
polarization states are taken into account. Note that long wavelengths cross the 



Hubble radius towards the end of the first infiation, so that H(tk) = ylm^ln^ in 
agreement with our numerical simulations 



3 32TrGml h 32']TGml r— 

k2h{k) = ]^ 3 V ^V^(^)-^0' ^>^o (10) 

where kb is the characteristic scale of the spectrum corresponding to the end of the 
first infiation driven by (ph, s = 2TiG{(f)f + 0^) is the number of e-folds till the end 
of the second infiation, s{k) = s{tk) = In = In fcj = ajHf where the index / 
denotes the values of all quantities at the end of the second infiation, and sq ^ 1 is 
the total number of e-folds during the second infiation. We take Sq = 27iG(f)l ^ 60, 
in order to get the characteristic scale kj, at the right place today. For the scalar 
adiabatic perturbations we have M 



3 6\/2nG / I \ \ 24:TiGml , 

^--m = + = ^ ^ V ^W(^W-^o) , ^ > ^0 . (11) 

As mentioned above, we are mostly interested in those values of — for which 
the intermediate stage, though it exists, is not well pronounced and deviations from 
monotonous behaviour of h{k) and $(fc) are concentrated in the region |s — SqI ~ 1- 
If the Hubble radius is crossed during the second infiation, we obtain the usual result 
for infiation driven by one scalar field 



,3,,,, 327iGmf r— ,3^,,, j2A7rGmf 

k-2h{k) = d --^^s{k), k-^^{k) = ^ s{k), s<So. (12) 
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This finally yields 



for both s > Sq and s < Sq, aside from the region |s — So| < ~ 1 around the 

point s = Sq. For s ~ sq, R{k) ^ 0.43. The behaviour of R{k) obtained numerically 
(see Figs.2-3) is in excellent agreement with (^). The slope of R{k) is: 

d\nR{k) 1, ^, 1 ^ 

-^.5K-„, + l) = — «1 (14) 



3 



where ns = I + = '^^"^ik.g^fc''^^'' • accordance with (14), we see 



from Fig.l that h{k) has a shape very similar to that of $(A;). However, this is not 



surprising since we have shown numerically [jT6| that the contribution of the light 
scalar field perturbations to $(fc) is dominant for the values of ^ considered here, 
while on the other hand the equation governing the light scalar field fluctuations 
essentially coincides with eq. (^, the equation for gravitational waves. 
We see that the relation 

riT ^ ^5 — 1 (15) 
(proposed, e.g., in [|l^) which has a very narrow region of applicability even for the 



case of single inflation, namely it is valid for power-law inflation only while it is not 
valid for chaotic inflation (see also the discusion in 0]), turns out to be valid for 
double inflation for large scales satisfying k <^ kb = kfe~^° (and not too small), 
as can be seen from (0) and also from the figures. As shown in [|12|, the best fit 
to observational data is given by 27ik^^ ^ (6 — 10)h~^Mpc, so the condition k <^ 
kb {s{k) > So) is certainly satisfied for scales probed by CMB anisotropy experiments 
on large and intermediate angles. On smaller scales for which k > kb, the relation ([T5|) 
is not valid any more. Hence, it cannot be considered as a test neither for generic 
inflation nor even for a single inflation. More general is the following relation which 
is valid for any slow-roll inflation driven by one scalar field: 



aw / o 
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ICi —R{k) = ICilnrl , < , (16) 



/C, = const = ^ (l + 1 ~ 6.20 , />1, In^l < 1 , (17) 
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/Ca ^ 6.93, /Cs ^ 5.45, /C5 ^ 5.10, /Cio ^ 5.30, Iut] < 1 (18) 

with /C5 being the minimal value for integer / < 40. In these expressions, only the 
Sachs- Wolfe effect is taken into account for both types of perturbations. Eq. (p!7| ) was 
first derived in [15], and (|l^) follows from results obtained in [^9[ (see also recent 



calculations in [^]). However, the quadrupole was not detected by COBE and it is 
usually excluded from the data; the asymptotic value (|1^ is not reached either for 
multipoles / = 3 — 15 measured by COBE. Taking into account the Doppler and the 
Silk effects for adiabatic perturbations yields a further decrease of /Cio by ~ 10%. 
Thus, it is better to use the value /C/ ~ 5 when estimating the GW contribution to 
the actual COBE data. 

The relation (^) is strongly violated in the case of double inflation for s{k) > sq 
because we then have: 

I = }Ci{0.3R{k)Y = ICi{s{k))'\ M = {s{k) - sor\ (19) 

therefore, T/S ICi\nT\. Note that relation (^) is not valid either for the R + E? 
inflationary model |l^, 0. So, (|1^) represents a decisive test for single inflation 



driven by a scalar field. Still this test is difficult because it first requires identification 
of the GW contribution to the ^ fiuctuations. Also, it cannot be used to predict the 
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amount of this contribution using the total COBE data since both sides of relation (|TB|) 
(in contrast to (|l^)) contain quantities referring to GW. 

For s ~ So, we get T/S = /CioSg^ ~ 0.08. As contributions to the ^ fluctuations 
from gravitational waves and adiabatic perturbations are statistically independent, 
the total rms value of ^ at large angles gets increased by a factor Jl + {T / S), i.e. by 



about 4%, due to the GW background. Therefore, in our model the ^ anisotropy on 
scales probed by COBE comes mainly from adiabatic perturbations. However, when 
comparing predictions of the model with observations the GW contribution has 
to be taken into account. 

Now we consider the dip in the R{k) curve (see Figs. 2-4) lying in the transition 
region between the two "plateaus" that arises for ^ > 15, i.e. when the intermediate 
period between two phases of inflation is sufficiently pronounced. Though, as can 
be seen from the figures, it constitutes a small effect even for ^ = 25, it is very 
interesting theoretically for two reasons. First, its very existence is due to the effect 
discovered in p]: the necessity of taking into account a contribution of the heavy 
scalar field quantum fiuctuations to the adiabatic perturbations in this range of scales. 
Indeed, in the absence of break between two phases of infiation (^ < 15), quantum 
fiuctuations of the light scalar field produce the dominant contribution to adiabatic 
perturbations over the region involved since 0z = 0o ^ 0/i ~ 0.5Mp. But if there 
exists a pronounced intermediate power- law stage {somfmj^'^ <^ 1), the contribution of 
the heavy scalar field fiuctuations becomes important in the transition region between 
the upper "plateau" and the lower one. It is given by the second term in the r.h.s. of 
eq.(4.18) of ref . [||: 



^ 571~ TV^^ '^o(y) sm (a-) + (27r(Ti) (-) In - j , 

/ 2\l/6 

cTi ^ 0.2875, a ^ 2.804, A;i ~ A; J I , ki ^ k ^ kb , (20) 

see eq. (4.32) in § for A; > kb- On the other hand (see [gT[ § and Fig.18.7, p.312 
in 



H), 



e^^hik) = ^6TTGmlln{j) (^^^ \ sin(a^)|, h ^ k ^ h . (21) 



The first, oscillating, term in the r.h.s. of (|20D arises from the light scalar field 
fiuctuations and is proportional to the r.h.s. of (pT]), while the second does not. In 
the absence of the heavy scalar field contribution, R{k) would be a smooth function 
of k. But the second term becomes important near the points where sin (c"^) = 0, or 
k = 1.12fci preventing $(/c) from going to zero in contrast to h{k). This explains why 
R{k) has a dip. For the values of ^ considered here, this effect is not big and can be 
calculated quantitatively with numerical methods only. Neverhteless, the appearance 
of the dip in the R{k) curve clearly shows that the heavy field contribution is present 
even if definitely subdominant. 
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A second new and remarkable point is that the minimum of the dip is not exactly 
zero as follows from the previously known eg . (^1]) . This is a direct consequence of 
the GW background being produced from the initial vacuum state during the first de 
Sitter stage via the process of particle creation, or in other words, of the GW being 
in a squeezed vacuum quantum state during inflation. The mathematical expression 
of this fact is that the amplitude hx{k) of the GW (operator) field modes has the 
following form in the regime k <^ aH: 



The first term in ( P2D is the quasi-isotropic mode which was made purely real by a 
time-independent phase rotation, the second one is the decaying mode. The condition 
Ci'^C2 = — I expressing quantum coherence between the two modes reflects the 
fact that initially, for t —>■ —oo, k ^ aH, the quantum state of the mode was the 
vacuum. Both terms in the right-hand side of eq. (^2]) are of the same order at 
the Hubble radius crossing during the first inflation. After that, the decaying mode 
quickly becomes exponentially small. If the latter mode is neglected as is usually 
done, then quantum coherence is lost and the resulting state of the GW for given 
k becomes equivalent to a classical wave with a Gaussian stochastic amplitude, a 
fixed temporal phase and some special correlation between k and — k modes (see 



the detailed discussion in ||23|). Actually, the quantum coherence is most probably 



lost shortly after the end of the second inflation due to interactions of the decaying 



mode with the environment though, as explained in there is no necessity to 
consider these interactions explicitly. As a result, the present GW background is not 
in a squeezed vacuum quantum state, it is purely classical. However, it is natural 
to assume that quantum coherence is not lost throughout inflation, at least if the 
decaying mode is not negligible. 

Then, for the double inflationary model, the decaying mode appears to be non- 
negligible over some range of k close to k^ due to the existence of a break between 
the two inflationary stages which yields two more Hubble radius crossings for modes 
belonging to that range. Keeping both terms in (|2^), in contrast to previous calcula- 



tions [^, ^, and using the method developed in [^], we obtain the following result for 
hx{k) after the last Hubble radius crossing during the second inflation (in the limit 
Somfmf^^ -C 1): 

k k 
hx{k) = Dism{—cr) + iD2COs(—a), 
ki ki 

= -5aV^, (23) 



3 H{tt) kf 

where Ho is the Hubble parameter at the beginning of the second inflation. Note that 
hx{k) could be made purely real had we chosen to make a different time- independent 
phase rotation. 

It is clear that the minimum of the dimensionless quantity k^^'^h{k), though small. 
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is not zero and it is reached at /c ~ vr 



1 . 



3 cr-^ H{tk) m,jMp 

hmaM 9 \k,) H\h) 9^6 ^2 1^(^1)-- ■ rnl ■ ^ ^ 

Hence, though the mode hx{k) is semiclassical after inflation, we may say that there 
appears a "quantum signature" in the present spectrum of the GW background due 
to the fact that the decaying mode should be taken into account in order to derive 
the correct form of the spectrum. 

Eq.(^) is in very good agreement with the numerical simulation displayed on 
Fig.l, where the dimensionless quantity k^^'^h{k) is expressed in uih/Mp units. For p = 
25, the value of the minimum ^ 0.2 which yields the value ~ 1.14, in remarkable 
agreement with (|2^) . Keeping the decaying mode in the scalar perturbations too, we 

get for Rrain- 

R • - ^° , ^ - ^0 . 1 (25) 

V'+V 8 ) V^+(, 8.3 ^2 ,3/2 

for the first minimum k = ki-. The value of ^r^^ and Rmin in the higher minima 
k = nki^, n > 1, quickly grows and the oscillations become unnoticeable. 

Thus, the spectrum h'^{k) has no zeroes. Note that the necessity of taking into 
account both modes resulting in the absence of zeroes (despite the existence of deep 
wells) in the perturbation spectrum may arise even for adiabatic perturbations gen- 
erated during inflation driven by one scalar field, provided the slow-rolling condition 
breaks for a short period of time, e.g. near values of the scalar field where the first 
derivative of its potential experiences a fast change So, we can formulate the 



following quantum positivity principle: if perturbations are generated by the process 
of particle creation from the vacuum, their power spectrum P{k) is non-zero for all 
modes. Of course, the reverse statement that the absence of zeroes implies that the 
fluctuations arise from quantum vacuum fluctuations, is not true. But the appearance 
of zeroes in a power spectrum implies the classical origin of perturbations, and the 
acoustic (Sakharov) oscillations in a purely baryon-radiation Universe is an example 
of a classical process leading to a power spectrum with zeroes for some values of k. 
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Figure 1: The quantities log k^^^^{k) (upper curve) and logk^^^h{k) (bottom curve), 
with the dimensionless quantities A;'^/^$(A:) and k^^'^h{k) expressed in ^ units, are 
displayed as functions of log A; (/cscale is arbitrary) for ^ = 25, k <^ aH. For the 

upper curve a{t) oc t^/^ nowadays is assumed. We note the presence of an oscillation 
in the transition region. The minimum of k^/'^h{k) is not zero due to the presence of 
the decaying mode as explained in the text. We note also the similarity of both curves 
for the /c-range displayed here corresponding to wavelengths that cross the Hubble 
radius towards the end of the first inflation (upper "plateau") and at the beginning 
of the second inflation (lower plateau) . 



Figure 2: The quantity logi?(A;) = log^^ is displayed as a function of log A; for 
^ — 25. The k scale is the same as in Fig.l. A dip is clearly seen in the transition 
region due to the contribution of the heavy scalar field fluctuations (see text) . Note 
again that the bottom of the dip is not zero. 



Figure 3: The quantity log R{k) = log is displayed as a function of log k for: 
q\ miL — 20 b) ^ = 15. Again, a dip is seen which decreases with decreasing 
due to the fact that for decreasing values of the intermediate stage becomes less 
pronounced and the contibution of the heavy scalar field fiuctuations decreases in the 
transition region. 



10 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9505125vl 



